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Abstract 

Simple, self-similar, analytic solutions of (l+3)-dimensional relativistic hydrody- 
namics are presented for ellipsoidally symmetric finite fireballs corresponding to 
non-central collisions of heavy ions at relativistic bombarding energies. The hydro- 
dynamical solutions are obtained for a new, general family of equations of state with 
the possibility of describing phase transitions. 
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1 Introduction 



Recently, new families of exact analytic solutions of relativistic hydrodynam- 
ics have been found in a Hungarian-Brazilian collaboration. The simplest case 
corresponds to the (l+l)-dimensional expansions [1]: it consists of a class of 
1-dimensional scaling flows with a proper-time dependent pressure, but char- 
acterized with an arbitrary rapidity distribution and a temperature field that 
is coupled to the rapidity distribution, thus overcoming one of the shortcom- 
ings of the well-known Hwa-Bjorken solution [2], which has a flat rapidity dis- 
tribution. Physically, the situation described in [1] may characterize e.g. soft 
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collisions of hadrons at high bombarding energies. These solutions are obtained 
for a broad class of equations of state that may include massive particles and 
an arbitrary constant of proportionality between the kinetic energy density 
and the pressure. The (l+l)-dimensional solutions have been generalized to 
the case of cylindrical symmetry in ref . [3] , describing a physical situation that 
may correspond to central collisions of heavy ions at ultra-relativistic collid- 
ing energies, thus, overcoming another shortcoming of Hwa-Bjorken solution, 
which contains no transverse flow. These exact analytic solutions, reported in 
refs. [1,3], are associated to each non-negative scaling function V(s) that sat- 
isfies the normalization condition V(0) = 1, where the argument s is a scaling 
variable that guarantees the self-similarity of the solutions. Here we present a 
broader class of solutions belonging to this new family of exact solutions. In 
particular, we shall describe two classes of ellipsoidally symmetric solutions of 
relativistic hydrodynamics, that may be relevant for the study of non-central 
collisions of heavy ions at relativistic bombarding energies. 

In general, exact analytic resolution of 3 dimensional relativistic hydrodynam- 
ics is a difficult task due to the highly non-linear nature of these equations. 
The complications are sometimes simplified by assuming idealized boundary 
conditions and simplified equations of state (e.g. infinite bombarding energy 
and a massless relativistic gas in the case of Hwa-Bjorken's solution [6]). A 
more realistic but analytically more complicated solution had been found by 
Khalatnikov [4], following Landau's basic ideas [5], that gave rise to the hy- 
drodynamical approach in high-energy physics. Both of these solutions are 
frequently utilized as the basis for estimations of various observables in ultra- 
relativistic nucleus-nucleus collisions [6]. 

Analytic solutions of relativistic hydrodynamics were also reported recently 
in refs. [7,8]. However, these solutions are valid only at the softest point of 
the equation of state. For more realistic situations, the (l+3)-dimensional 
relativistic hydrodynamical equations are frequently solved using various nu- 
merical methods, for example, recent solutions were obtained with the help of 
smoothed-particle hydrodynamics in refs. [9,10]. 

The exact analytic solutions, reported here, are generalizations of the results of 
refs. [1,3] to the case of 3-dimensional relativistic expansions with less and less 
symmetry. However, we emphasize that these results also correspond to gener- 
alizations of earlier analytic solutions of non-relativistic hydrodynamics. The 
first solutions of this kind have been found by Zimanyi, Bondorf and Garpman 
(ZBG) in 1978 [11]. The key aspect of the ZBG solution is its self-similarity, 
this property is kept in all the subsequent generalizations including cylin- 
drical or ellipsoidal symmetries and relativistic flows. The ZBG solution has 
been generalized to the case of spheroidally symmetric expansions in ref. [12], 
to spherically symmetric Gaussian expansions in ref. [13]. From the point of 
view of the relativistic generalizations, an important stepping stone was made 
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in ref. [14] when a scaling function has been introduced for the spherically 
symmetric non-relativistic solutions of hydrodynamics. This was the first case 
when a whole new family of solutions has been found, that assigned an inde- 
pendent solution for every (differentiable and integrable) scaling function that 
satisfied the requirement of positivity T(s) > and normalized as T(0) = 1. 
The ZBG and the spherical Gaussian solutions appeared as special choices 
for the functional form of the scaling function T(s). The Gaussian family of 
spherical solutions has been generalized to Gaussian ellipsoidal expansions in 
refs. [15,16]. Recently, the whole family of self-similar ellipsoidally symmetric 
non-relativisitic solutions of hydrodynamics has been found in ref. [17]. In the 
present paper, we report on the generalization of this family of solutions to 
the case of relativistic expansions with ellipsoidal symmetry and introduce 
new families of relativistic solutions with even less than ellipsoidal symme- 
try. All the solutions reported here have non-trivial non-relativistic limiting 
behaviour. 



2 The equations of relativistic hydrodynamics. 



Let us adopt the following notational conventions: the coordinates are x^ = 
(t, r) = (t,r x ,r y ,r z ), x^ = (t, —r x , —r y , —r z ) and the metric tensor is g^ v = 
g^ u = diag(l, —1,-1,-1). The relativistic continuity and energy-momentum- 
conservation equations are 



d>o = o, (i) 

d u T" u = , (2) 

where n = n(x) is the conserved number density, w M = u^(x) = 7(1, v) is the 
four- velocity, with u^u^ = 1, and T^ v is the energy- momentum tensor. We 
assume perfect fluid, 

= (e + p)u^u v - pg^ , (3) 
where e stands for the relativistic energy density and p is the pressure. 

From general thermodynamical considerations and the perfectness of the fluid 
one can show that the expansion is adiabatic, 

9>^) = 0, (4) 

where a stands for the entrophy density. 

For the equations of state, we assume a gas that may contain massive conserved 
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quanta, 



e = mil + Kp , (5) 
p = nT. (6) 



These equations have two free parameters, m and k. Non-relativistic hydro- 
dynamics of ideal gases corresponds to the limiting case of v 2 <C 1 and k — | . 
Our solutions, presented below, exist for all values of m > and k > . 



Using the continuity equation (1) and the equations of state (5,6), the energy- 
momentum-conservation equations, (2) can be transformed to the Euler and 
temperature equations, 

u^dup + (e + p)u»d^u v - d v p = , (7) 
u^df/T + W = . (8) 



Expressing the energy density and the pressure in the 3 independent compo- 
nents of the relativistic Euler equation in terms of n and T with the help of 
the equations of state, (5,6), one obtains a closed system of 5 equations (the 
continuity, the Euler and the temperature equations) in terms of 5 variables, 
n, T and v = (v x ,v y ,v z ). 



3 Self-similar, ellipsoidally symmetric ansatz 



We search for self-similar solutions in which the isotherms at each instant 
are ellipsoidal surfaces, where the number density is also constant. We do 
not discuss here how the major axes of these ellipsoids may be rotated in 
the frame of observation. We describe the solution in the natural System of 
Ellipsoidal Expansion (SEE) [15,16], where the coordinate axes point to the 
pricipal directions of the expansion. Let us now define the following scaling 
variable 

s = — + + — (9) 

X 2 Y 2 Z 2 ' 

where the scale parameters X — X(t),Y — Y(t) and Z = Z(t) are assumed to 
depend only on the time variable t. The condition s = const., at each instant 
t, defines the isotherms mentioned above. As these surfaces are ellipsoids, the 
solutions belonging to this class will be characterized by ellipsoidal symmetry. 
We assume a self-similar expansion of Hubble type, with possibly different 
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Hubble constants in all the pricipal directions of the expansion: 



v x (t,r) = — r x , (10) 
Y 

v y(t, r ) = y Ty ' ^ 1V) 

v z (t,r) = ^r z . (12) 

where A = dA(t)/dt stands for the derivative of the scale parameter A = 
{X, Y, Z} with respect to time t. 

As a straightforward generalization of the results derived in refs [1,3], we find 
the following new family of ellipsoidally symmetric, exact analytic solutions 
of relativistic hydrodynamics: 

™2 ^2 2 

' +7^7 + ^, (13) 



X 2 t 2 Y H 2 Z 2 t 2 



V 


r 

~~ t 




ra(t,r) 


= n | 


Vr 


p(t,r) 


= Po ( 


'TO 

\ T , 


T(t,v) 


= T \ 


fro 



3+3/k 



or u» = —, (14) 

T 

(15) 
(16) 
(17) 



v( s y 

where r = y/t 2 — r 2 , p = n T , X , Y , Z are constants and V(s) is an 
arbitrary positive, differentiable function of the scaling variable s normalized 
such that V(0) = 1. The temperature distribution also depends on the scaling 
variable s through a scaling function T(s) , which happens to be T(s) = 
1/V(s) . Note that the pressure p depends only on r. 

An interesting property of this solution is that the flow and the pressure 
fields are spherically symmetric, however, the other thermodynamical quan- 
tities such as the temperature or density distributions are characterized by 
ellipsoidal symmetry. 

The lack of acceleration is reflected by the equations 

X = X t, (18) 
Y = Y t, (19) 
Z = Z t . (20) 

In this family of new solutions, the flow is three-dimensional, accelerationless 
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Hubble-type, = x^/r, which gives 



d^u v = 0. 



(21) 



This property, together with the s-independence of pressure fields, guarantees 
that the Euler equation be satisfied regardless of the mass m and the value of 
k. We have found new, self-similar, scaling solutions to the continuity and the 
temperature equations with ellipsoidal symmetry. Thus a new hydrodynamical 
solution is assigned to each scaling function V(s) = 1/T(s), similarly to the 
cases of the ellipsoidally symmetric, non-relativistic solutions of ref. [17]. 



4 Phase transitions and freeze-out 

The simplest case of the previously discussed new family of exact solutions 
of relativistic hydrodynamics is given by the choice of V(s) = 1. In this case, 
the temperature is constant on a hypersurface characterized by a constant 
proper-time, and all the other thermodynamical parameters are constants on 
these surfaces as well. Thus these parameters depend on r only, similarly to 
the Hwa-Bjorken solution [2]. 

It is possible to generalize the equations of state, (5,6)to characterize a new 
phase of matter and then describe the rehadronization of a "Quark Matter" 
(QM) state to a Hadron gas (H) in a model similar to that of Gyulassy and 
Matsui [18]. For simplicity, from now on we refer to the the equations of state 
(5,6), describing massive, conserved quanta, as that of a hadron gas and index 
the variables with subscript h- To phenomenologically describe a new phase, 
which includes constituent quarks and anti-quarks (Q), characterized by their 
mass and by a vacuum pressure, we generalize the equations of state to 



where \ £ and A p are constants. Let us introduce the notation kq = \ £ /\ p . 
With this simple ansatz, the form of the relativistic Euler and temperature 
equations does not change and we obtain the same form of the solutions for n, 
T and v as before, however, with modified boundary conditions corresponding 
to modified constants of integration. We may refer to this new phase as a kind 
of quark matter, providing simple model equations of state for massive quarks, 
and a bag constant corresponding to color deconfinement. The gluons are 
assumed to be integrated out, providing mass for the quarks. Such a picture is 
qualitatively supported by phenomenological fits to the lattice QCD equation 
of state with quarks that pick up the constitutent mass around T c and gluons 
that carry big effective mass in the same temperature domain [19]. 



E Q = vtlqUq + X £ n Q T + B, 
Pq = X p n Q T - B, 



(22) 
(23) 
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Using the above ansatz for the equations of state, we obtain the following 
solution for the pressure: 



If we chose B > and \ p > 1, then at low temperatures the state consist- 
ing of non-relativistic massive ideal gas will be stable, and at some critical 
temperature T c the pressure of the two phases is balanced. However, at high 
temperatures, the pressure of the new phase will be larger than that of the 
hadronic matter, as \ p > 1, hence at high temperatures this phase will be the 
stable one. 

The critical temperature, T c is defined by 



For the hydrodynamical description, the specification of Tin and Uq as a func- 
tion of e.g. a temperature and chemical potential is not necessary, any ther- 
modinamically consistent form can be included into the hydro solutions that 
we focus on. Hence the above equation may indicate, if nn and hq are ex- 
pressed as a function of the temperature and the (baryo) chemical potential, 
that the critical temperature may become density dependent. 

One can show that the expansion is adiabatic if the matter consists of only one 
of the phases. We may assume, that the expansion is adiabatic also during the 
period of phase coexistence, and describe the phase transition with the help 
of a Maxwell construction. This case is similar to the construction discussed 
by Gyulassy and Matsui in case of a massless ideal gas [18]. Let us denote the 
entrophy density of the two phases by Sq and sh, and the volume fraction 
of the two phases by Jq and fn = 1 — /q, respectively. The proper-time 
dependence of the entrophy density is then given by 



At the beginning of the phase transition, at r = tq c the whole system consists 
of "Quark Matter", that is /q(tq iC ) = 1. By the time the phase transition ends, 
one obtains /q(t# iC ) = 0. As the well known steps of the Maxwell construction 
can be copied from the case without conserved charges, for further details we 
simply refer to section 6.2.3 of ref. [6]. 

Starting the hydrodynamical expansion from a quark matter initial state, the 
form of our hydrodynamical solution becomes unmodified even during the 
phase transition, and the coasting type of hydrodynamical evolution can be 
continued at a fixed T c , from tq iC to th, c untill all matter is transformed to 
hadronic matter. Then one can continue the same kind of coasting solution 
untill the hadronic matter freezes out. 




(24) 



n H T c = X p n Q T c - B. 



(25) 



s ( r ) = /q(t)sq,c + (1 - fQ(r))s H , c - 



(26) 
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This expansion takes usually a long time, which is difficult to justify by ex- 
perimental evidences. However, the same kind of hydrodynamical solution 
can be matched with a non-equilibrium type of phase transition, the sudden 
rehadronization and a simultaneous freeze-out of hadrons from a supercooled 
Quark Matter, as described in greater details in refs. [20,21]. Then in the same 
model one may reach a deeply supercooled Quark Matter state, and the tran- 
sition to the hadron gas may proceed via a mechanical instability, associated 
with a negative pressure state, governed by the conservation of matter, energy 
and momentum and the impossibility of an entrophy decrease in a deflagration 
through a hypersurface with a time-like normal vector at r = r TD , 



where dfn v = da v is the normal- vector of the r = ttd hypersurface. These 
equations and inequality can be solved very similarly to ref. [20] , that discussed 
the case of a transition from a massless ideal quark-gluon plasma (QGP) to 
a massless hadron gas. The lack of acceleration, and the invariance of the 
equations of state for rescaling the temperature and for adding a vacuum 
energy term are the essential reasons why earlier considerations that discussed 
phase transitions for massless particles can be straigthforwardly implemented 
to our new family of exact hydrodynamical solutions that describe the phase 
transitions and expansions of massive quanta. 

Recently, a new method has been proposed to evaluate particle emission from 
hydrodynamically evolving, locally thermalized sources, that generalizes the 
Cooper- Frye freeze-out conditions for systems with volume emission [23,22]. 
This method, based on the method of escaping probabilities, can also be ap- 
plied to evaluate the particle spectra and correlation functions from the new 
hydrodynamical solutions, as the non-relativistic form of the solutions corre- 
sponds to the limiting behaviour that was discussed in ref. [23]. 



5 Factorized solutions 



From the structure of the new family of solutions presented above, it is obvious 
how to generate further new solutions of relativistic hydrodynamics. 

Observe that a modified version of the scaling variable of eq. (13) can be 
defined, and another family of the same kind of ellipsoidal solutions can be 
generated by introducing 



[T^dfn v ] = 0, 
[nu u dfn u ] = 0, 
[su u dfn u ] > 0, 



(27) 




(29) 
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,y*2 ^2 ^2 

s ' = xp + Yp + zp' (30) 

X' = X Q r, (31) 
V = Y t, (32) 
= Z t, (33) 

where the dots should be understood as derivatives with respect to r. Mathe- 
matically, the modified scaling variable solves the relativistic hydro equations, 
because it also satisfies the requirement 



u^s = u^d^s' = , (34) 

if the flow is scaling. This condition can be considered also as the criterium of 
a "good" scaling variable, and satisfied by the following independent sets of 
variables: 

s x = T f, (35) 
Sy = -f , (36) 



r 2 

2 



z t 2 ' 

or alternatively, 



(37) 



4 = 4, W) 



T 2 



rl 



4 = ^ (39) 



4 = r f 2 ■ (40) 



All of the 6 scale variables defined above satisfy eq. (34) if the flow profile is 
of a 3-dimensional scaling type, = x^/r . 

As the relationship V(s) = 1/T(s) holds in the class of solutions we present 
here, if we assume a factorized form for the scaling function of the density, 
we automatically generate a factorized form for the scaling function of the 
temperature, hence it becomes easy to generate further new solutions. Let us 
define the scaling functions V x (s x ), V y (s y ) and V z (s z ) that are positive, differ- 
entiate and satisfying 14(0) = V y (0) — V z (0) = 1, otherwise being abritrary 
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and independent from each other. Then, a new type of hydro solutions reads 



V 


t 




n(t,r) 




It 


(t,r) 


= Po [ 


'TO 
V T , 


T(t,r) 


= T \ 
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or v} 1 = —— , (41) 



V x (s x )V y (s y )V z (s z ) , (42) 

3+3/k 

x3/k 1 I I 

r) V x (s x )V y (s y )V z (s z ) • (44) 



(43) 



Note that this form of solution is invariant for a change of the scaling variables 
and the scales as (s x ,s y ,s z ) — > (s^s^s^,) an d (X,Y,Z) — > (X', F' , Z'), with 
the time derivatives in the definitions of the Hubble flow field, eqs. (10-12) 
understood as derivations with respect to r as in eqs. (30-33). 



6 General form of solutions, beyond ellipsoidal symmetry 

The key point in checking that the above factorized forms indeed solve the 
equations of relativistic hydrodynamics is that the comoving derivatives of the 
scaling functions, u^d^V^Si), are proportional to u^d^Si , hence vanish. This 
is, together with u^d^Uy = 0, the essential property of the scaling solutions. 

We find that the most general form of the scaling variable is 

s = F(s x ,s y ,s z ) = G(s' x ,s' y ,s' z ), (45) 

which means that any function of scaling variables (s x , s y , s z ) or (s' x , s' y , s' z ) 
can be utilized as a new scaling variable. Indeed, we have 

BF 

u^F(s x ,s y ,s z )= ]T —u^ Si = Q. (46) 

i=x,y,z l 

which yields the generalized form of the new family of solutions of relativistic 
hydrodynamics: 





r 


V 


~ t 


n(t,r) 


= n 


p(t,r) 


= Po 


T(t,r) 


= T 



T 



or u» = — , (47) 

r 

V(s) , (48) 

_ x 3+3/k 

7) • 
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with the constraint that V(0) = 1. For example, if we choose s = s x + s y + 
s z , we obtain solutions with ellipsoidal symmetry, however, one may choose 
s — s x + s y — s z to obtain solutions where the isotherms are one-sheeted 
hyperboloids, or s = —s x — s y + s z where the temperature and the density are 
constants on two-sheeted hyperboloids. In fact the possibilities are infinitely 
rich. 



7 Non-relativistic limiting behaviour 



As the considered equation of state contains mass as a free parameter, it is 
possible to study the non-relativistic limiting case of these new relativistic 
solutions. The local thermal motion is non-relativistic if m > T, the flow 
is non-relativistic in the region of |r| <C t, which implies r ~ t and relates 
the relativistic solutions with ellipsoidal symmetry to the asymptotic, large 
time behaviour of the ellipsoidally symmetric, non-relativistic solutions of 
refs. [15,16,17]. Let us recapitulate here the general form of the non-relativistic 
solutions of fireball hydrodynamics [17]. 

The definition of the scaling variable s coincides with that of eq. (9). The 
solution for the flow field coincides with that of eqs. (10-12). The scaling 
function for the temperature, T(s) is utilized to express the solution for the 
NR density and temperature fields as 



n 



n - 



XqYqZq 1 
XYZ T(s) 
XqYqZq 
XYZ 



cxp 



m(X 



+ Y 2 + 7 2 ) 

as ' as ' as) 



T 



du—— 
T(u) 



1/k 



T{s). 



(51) 
(52) 



The large time behaviour of these non-relativistic, self-similar, ellipsoidal so- 
lutions of hydrodynamics is characterized by constant values of (X as , F as , Z^), 
which implies that for asymptotically long times a scaling flow field develops 
also in these non-relativistic solutions, v ~ r/t. In the exponential factor that 
appears in this non-relativisitic solution for n we may expand the inverse scal- 
ing function 1/T(u) as a polinomial in u. Keeping the leading order terms, 
and performing the integration, we find that the exponential factor will be 
approximately a Gaussian factor. However, the widths of the Gaussians in 
all the directions will be proportional to the time. As we have assumed that 
|r| <^ t, this exponential factor yields a factor of l + (9(r 2 /t 2 ), hence we find the 
general asymptotic behaviour of the non-relativistic solutions to be the same 
form as the general solution of the relativistic solution of the hydrodynamical 
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equations: 



n, 



as 




(54) 



(53) 



Due to the rescaling, the constants of normalization n as and T as are different 
from n and T . 



8 Summary 

We have generalized the recently found new family of solutions of relativistic 
hydrodynamics to the case of expanding fireballs with ellipsoidal symmetry. 
The solutions contain an arbitrary scaling function V(s), restricted only by 
non- negativity and by the requirement of V(0) = 1, a very rich set of possible 
scaling variabes s, and 5 important parameters, the mass m, the parameter 
k of the equation of state, the scale parameters X , Y and Z . Furthermore, 
we generalized the equations of state to describe a phase transition from a 
deconfined state consisting of massive (constituent) quarks and antiquarks to 
a state consisting of massive hadrons. Based on the scaling properties of the 
equations of state, as well as on the coasting nature of the expansion, we 
have shown that both the usual adiabatic Maxwell construction as well as the 
fast, non-equilibrium time-like deflagrations can be constructed and described 
within the considered class of hydrodynamical solutions. 

We expect that some of these solutions of relativistic hydrodynamics may 
have future applications in non-central collisions of heavy ions at relativistic 
bombarding energies, a topic of great current research interest. 
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